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ABSTRACT 

The effect of the strange quark in the finite temperature phase transition of QCD 
is studied on the lattice. Using the one-plaquette gauge action and the Wilson 
quark action, the transition in the chiral limit is shown to be continuous for the 
case of degenerate two flavors, Np = 2, while it is of first order for Np > 3. 
For a more realistic case of massless up and down quarks and a light strange 
quark, Np = 2 + 1, clear two state signals are observed both for m s ~ 150 and 
400 MeV. In contrast to a previous result with staggered quarks, this suggests 
a first order transition in the real world. In order to see the implication of these 
results to the continuum limit, we started to study these issues using improved 
actions. First results using a RG improved gauge action combined with the 
standard Wilson quark is presented for the case of Np = 2: With this action 
the finite temperature transition is shown to be continuous in the chiral limit 
confirming the result of the standard action. Furthermore, not like the case of 
the standard action where lattice artifacts make the transition once very strong 
at intermediate values of the hopping parameter K on N t = 4 and 6 lattices, a 
smooth crossover is found for the improved action when we increase 1/K — 1/ K c , 
in accord with a naive expectation about the fate of second order chiral transition 
at finite m q . 

1. Introduction 

At sufficiently high temperature, QCD looses two of its characteristic properties: 
the confinement and the spontaneous breakdown of chiral symmetry. Studies on the 
lattice have shown that the nature of the corresponding phase transition/crossover 
sensitively depends on the number of flavors Np and quark masses: With degenerate 
Np quarks, the transition is first order in the chiral limit for Np > 3, while a contin- 
uous transition is found for Np = 2. At sufficiently large quark mass the transition 
is smoothed away into a crossover. Therefore the nature of the transition in the real 
world may depend sensitively on the values of m ? 's, especially on the s-quark mass m, 
because m s is of the same order of magnitude as the estimated transition temperature 
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T - 100 - 200 MeV, 

In this paper, I report our recent study of finite temperature QCD transition 
on the lattice. We use Wilson fermions as lattice quarks. The investigation using 
Wilson quarks is highly important because the Wilson formalism of fermions is the 
only known formalism which possesses a local action for any number of flavors with 
arbitrary quark masses. We mainly perform simulations on lattices with the temporal 
extension N t = 4. Although N t = 4 is not large enough to obtain the result for the 
continuum limit, this work is a first step toward the understanding of the nature of 
the QCD phase transition with Wilson quarks. 

We first study this using the standard one-plaquette action for the gauge part. 
The iVi?-dependence of the chiral transition is discussed in the next section. 1 I then 
discuss the influence of the s-quark for the case of degenerate three quarks (Np = 3) 
and the case of two massless u and d-qaurks with a light s-quark (Np = 2 + l). 2 

In order to see the implication to the continuum limit, a careful estimation of 
lattice artifacts is required. This motivated us to study these issues using improved 
actions. In Sect. 3, first results using a RG improved gauge action combined with 
standard Wilson quark is presented for the case of Np = 2. 3 ' 4 

2. Chiral transition 1 

In the coupling parameter space (f3,K) of lattice QCD, where (3 = 6/g 2 is the 
gauge coupling and K is the hopping parameter, the location of the chiral limit K c 
is defined as the point where pion mass m v vanishes in the chirally broken phase. 
Alternatively we can define K c as a point where the quark mass m q vanishes, 5 ' 6 
with m q defined through an axial vector Ward identity. 7 ' 8 The location of the finite 
temperature transition/crossover K t can be easily identified numerically by rapid 
changes of physical observables such as the plaquette, the Polyakov loop and hadron 
masses. Our study has shown that, for Np < 6, the line K t crosses the line K c at 
finite (3. 6,1 Phase diagram for Np = 2 and 3 is given in Fig. 1. 

In order to study the nature of the transition at the crossing point f3 ct of K c and 
K t , we perform simulations on the line K c ((3) starting from a (3 > (3 ct and reduce (3. 
The number of iterations A^; nv needed for the quark matrix inversion provides a good 
indicator to discriminate the deconfining phase from the confining phase, especially 
on K c , because 7Vj nv is enormously large on the K c line in the confining phase while 
it is of order several hundreds in the deconfining phase. This difference is due to the 
appearance of zero modes around K c in the confining phase 9 . Combining this with 
measurements of physical observables, we identify the crossing point (3 ct . We call this 
method as "on-A^" simulation method. 1 ' 4 We can check that f3 ct thus determined is 
consistent with an extrapolation of the K t line toward the chiral limit (Fig. 1). From 
the behavior of physical quantities toward p ct , we are able to determine the order of 
the chiral transition. 
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Fig. 1. Phase diagram for Np = 2 and 3. Filled symbols are for K c (rn^) and K c (ra q ). Open 
symbols are for K t (N t = 4) for Np = 2 while open circles with cross for Np = 3. Recent results 
of various groups are collected (see references in 4). Circles (filled, open and those with +) are our 
data. Shaded lines are to guide the eye. 
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Fig. 2. (a) Time history of 7V; nv and (b) /3-dependence of (m^a) 2 for Np = 2 on K c obtained on an 
8 2 x 10 x 4 lattice. 
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Fig. 3. Time history of iV; nv for Np = 3 on K c 
obtained on an 8 2 x 10 x 4 lattice. 
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Fig. 4. (m^a) 2 for m s ~ 0, 150 and 400 MeV 
with m u d — 0. Filled and open symbols are for 
8 2 x 10 x 4 and 12 3 x 4 lattice, respectively. 



Time history of 7Vj nv for Np = 2 on the ii^-line is shown in Fig. 2(a) for N t = 4. 
We find that when (3 > 4.0, 7Vi nv stays around several hundreds, while (3 < 3.9 it 
increases rapidly with r. Time histories of the plaquette, the Polyakov loop, m v etc. 
confirm that the system is developing into a confined state for (3 < 3.9. Therefore 
we conclude that f3 ct ~ 3.9 - 4.0 for N t = 4. Similar study on an N t = 6 lattice 
gives ~ 4.0 - 4.2. For both N t , we find no two-state signals around (3 ct . m% on 
the i^ c -line is shown in Fig. 2(b). We find that, when we decrease f3 toward ct , 
decreases to zero smoothly. This suggests that the transition for Np = 2 is continuous 
in the chiral limit, in accord with the result of an universality argument 10 and scaling 
studies with staggered quarks. 1112 

Corresponding results for Np = 3 are given in Fig. 3 and Fig. 4. (In Fig. 4raj~0 
corresponds to the case of iV^ = 3 on K c . Results for Np = 2 + 1 will be discussed 
in the next section.) In sharp contrast with the case of Np = 2, we here find a clear 
two-state signal at (3 = 3.0 ~ (3 ct : For a hot start, 7Vi nv is quite stable around ~ 800 
and m 2 is large (~ 1.0 in lattice units). On the other hand, for a mix start, A^ nv 
shows a rapid increase with r and exceeds 2,000 in t ~ 20, and in accord with this, 
rr? v the plaquette, the Polyakov loop etc. decreases with r. This implies a first order 
chiral transition for Np = 3. These results are again consistent with the prediction 
based on universality 10 . 

3. The effect of the strange quark 2 

As discussed in the introduction the mass of the s-quark may play a decisive role 
on the nature of the QCD transition. Let us begin with the case of degenerate Np = 3 
and increase m q by increasing (3 along K t shown in Fig. 1. On the iQ-line we observe 
clear two-state signals for (3 < 4.7 (Fig. 5), while no such signals for (3 > 5.0. At 
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5. Time history of the plaquette for 
= 3 at f3 = 4.7 obtained on a 12 3 x 4 
lattice. 
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Fig. 6. Time history of the plaquette for 
N F = 2 + 1 at m ud ~ and m a ~ 400 MeV 
obtained on a 12 3 x 4 lattice. 



the transition point (in the confining phase) of ft = 4.7 we get m q a = 0.175(2) and 
(m v /m p ) = 0.873(6). Using our estimation of a -1 ~ 0.8 GeV for these ft, the critical 
quark mass m"' 1 up to which the first order phase transition persists can therefore 
be bounded from below as m c I lt > 140 MeV. a 

We note that these values are much larger than that obtained previously with 
staggered quarks 13 ' 14 : m^a = 0.025 - 0.075. Using a" 1 ~ 0.5 GeV at ft ~ 5.2 
obtained by an extrapolation of the data for N F = 2 15 we find m c J' lt ~ 10 - 40 MeV.' J 
Using a result of hadron masses at these ft for N F = 4 16 (b ecause the data for Np = 3 
is not available), these m r q nt correspond to (m v /m p ) crit ~ 0.42 - 0.58. 

We now study the case of Np = 2 + 1: m u = raj < m s . When m s is reduced 
from oo to 0, the phase transition must change from continuous to first order at some 
quark mass m*. In this connection, we note that, when m s < m^" lt (the critical 
quark mass for degenerate Np = 3), we expect a first-order transition, while when 
m s > m u< i > m^"'*, we expect a crossover. 

We study two cases m s ~ 150 MeV and 400 MeV. We apply the "on-if c " simula- 
tion method: Keeping K U( i = K c , i.e. m u d — 0, and K s at the point corresponding to 
the cases m, ~ 150 or 400 MeV, we decrease ft until we hit K t . Fig. 6 is an example 
of time history we get. We find two state signals for both m s = 150 and 400 MeV. 2 

on the i^-line is shown in Fig. 4. 

These results imply that m* > 400 MeV. With staggered quarks Columbia group 14 
reported that no transition occurs at m^a = 0.025 and m s a = 0.1 (m ut i ~ 12 MeV, 
m s ~ 50 MeV using a -1 ~ 0.5 GeV). Albeit slightly heavier masses for u and d 

a Our results of the hadron spectrum in the range of (3 = 3.0 - 4.7 for Np = 2 and 3 show that 
physical s-quark mass determined from = 1020 MeV is about 150 MeV for our definition of m q . 
h a~ 1 of staggered quarks for Np = 3 may be slightly (about 0.1 GeV) larger. This small difference 
does not affect our discussions below. 
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Fig. 7. Order of the finite temperature QCD transition in the (m u d,m a ) plane. First order signals 
are observed at the points marked with filled circle, while no clear two state signals are found at the 
points with open circle. The values of quark mass in physical units are computed using a -1 ~ 0.8 
GeV. The real world corresponds to the point marked with star. The second order transition line is 
suggested to deviate from the vertical axis as m u( i oc (m* — m s ) 5 / 2 below m*. 

quarks, this suggests a discrepancy between two lattice quarks. Our results on the 
nature of the transition for degenerate Np = 3 and Np = 2 + 1 are summarized in 
Fig. 7. 

4. Improved action 

In the previous sections we have seen that, although the results of Wilson quarks 
and staggered quarks are qualitatively consistent, quantitative predictions about the 
values of critical quark masses are discrepant yet on N t = 4 lattices. We need to get 
closer to the continuum limit to resolve the discrepancy. This motivated us to study 
these issues using improved actions. 

Based on a study of a block transformation, Iwasaki proposed a RG improved 
gauge action about ten years ago 17 

S g = l/g 2 {co y^(plaquette) + c\ y^(l x 2 loop)} 

with c\ = —0.331 and cq = 1 — 8ci. c As a first step toward an improved full QCD, 

c This choice of coupling parameters minimizes the distance to the renormalized trajectory after one 
block transformation. A perturbative calculation shows that the scale parameter of this improved 
action Ajm is close to that in the MS scheme: A-j^/A/m — 0.488, 18 in contrast with the large 
ratio for the case of the standard action. This action was applied to numerical studies of the string 
tension, the hadron spectrum, topological properties and the U(l) problem 19 . 




Fig. 8. (ro,fl) 2 from an improved Np = 2 
action on K c obtained on an 8 3 x 4 lattice. K c 
is determined by the vanishing point of ra\ on 
an 8 4 lattice. 
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Fig. 9. (m^a) 2 from an improved Np = 2 
action as a function of 1/K — 1/K C obtained 
on an 8 3 x 4 lattice. 



we combine this RG improved pure gauge action with the standard Wilson action. 3,4 
This should work well at small K where the effects of quarks can be absorbed by 
renormalizations of pure gauge interactions. Our numerical results suggests, however, 
that our action improves the theory very much also near K c . 3,4 

m% on K c is shown in Fig. 8. We see that m 2 decreases smoothly to zero when 
we decrease (3 toward (3 ct ~ 1.4 - 1.5. Therefore, as in the case of the standard 
action, the finite temperature transition is continuous in the chiral limit. Shown in 
Fig. 9 is the 1/K — 1/K C (~ m q ) dependence of m 2 for various (3. The straight line 
envelop in the right corresponds to the chiral behavior oc m q in the confining 
phase and deviation from this line signals the transition to the deconfining phase. We 
find that m 2 , as well as other physical observables, is quite smooth at large p. This 
strongly suggests that the transition there is actually a crossover, in accord with a 
naive expectation that the second order chiral transition will become a crossover at 
finite m q . This result is in clear contrast with that with the standard action where 
lattice artifacts make the transition once very strong at intermediate K on N t = 4 
and 6 lattices. 20 



5. Conclusions 



We have studied the nature of the finite temperature transition in QCD with 
Wilson quarks for Np = 2, 3 and 2+1. Using the standard one-plaquette action for 
the gauge part, the chiral transition is shown to be continuous for Np = 2 on both 
N t = 4 and N t = 6 lattices. For Np = 3 at N t = 4, clear two state signals are 
observed on the K f -\me for m n < 140 MeV. For Np = 2 + 1 we have studied the cases 
m s ~ 150 and 400 MeV with m u = rrid — 0, and we have found two state signals for 



both cases. These results suggest a first order QCD transition in the real world. Our 
results on the nature of the transition for degenerate Np = 3 and Np = 2 + 1 are 
summarized in Fig. 7. 

We also studied these issues using a RG improved action. As in the case of the 
standard action, the finite temperature transition is continuous in the chiral limit for 
Np = 2. On the other hand, not like the case of the standard action where lattice 
artifacts make the transition for Np = 2 very strong at intermediate K, the transition 
with the improved action becomes quickly smooth when we increase 1 /K—l / K c . This 
suggests that the transition there is a crossover, in accord with a naive expectation 
about the fate of second order chiral transition at finite m q . Studies for Np = 3 and 
2 + 1 are in progress. 
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